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Abstract - We compute the decay of the autocorrelation function of the observable 
\v x | in the Sinai billiard and of the observable v x in the associated Lorentz gas with 
an approximation due to Baladi, Eckmann and Ruelle. We consider the standard 
configuration where the disks is centered inside a unit square. The asymptotic decay 
is found to be C(t) ~ c(R)/t. An explicit expression is given for the prefactor c(R) as 
a function of the radius of the scatterer. For the small scatterer case we also present 
expressions for the preasymptotic regime. Our findings are supported by numerical 
computations. 



1 Introduction 

In this paper we are going to present analytical and numerical evidence that typical correlation 
functions in a Sinai billiard with infinite horizon asymptotically decay as C(t) ~ c(R)/t. On the 
theoretical side there has been some heuristic argument ^| and rigorous bounds (^] suggesting 
this type of behavior. We will strengthen the argument for the suggested type of decay and derive 
an expression for the constant c(i?) as a function of disk radius R. On the numerical side we 
investigate correlation functions via Monte Carlo simulations: a remarkably good agreement is 
found with theoretical expressions, which accurately predict the observed long-time behaviour, 
and also account for the main features of the preasymptotic region. The 1/t decay of the velocity 
autocorrelation function was, to our knowledge, first observed numerically in jy], and confirmed, 
in a somehow different approach, focusing on power spectrum behaviour, in 0. More recently 
the asymptotic algebraic decay was detected in |Q, while only the preasymptotic regime was 
investigated in the simulations described in 0. 

We will focus on two particular correlation functions. First the autocorrelation function of 
the observable \v x \ in the Sinai billiard. The modulus is taken to prevent the observable to 
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change sign by the frequent bounces on the square walls which would make our calculations more 
difficult. Secondly we study the autocorrelation of the observable v x in the associated Lorentz 
gas, the difference from the former case is that the sign of the observable may now be changed 
by bounces on the disk(s). 

The theoretical calculations will be carried out following an idea originally due to Baladi, 
Eckmann and Ruclle |j, [?], |[ g, [To). This approximation (called the BER approximation) for 
intermittent systems divide the temporal evolution of the systems into intervals, each interval 
consist of one laminar and one chaotic interval. The approximation assumes the absence of 
correlations between disjoint intervals. In particular the lengths Aj are assumed mutually in- 
dependend and one can define a probability distribtuion p(A), a function that encodes a lot of 
information about the system. For the observables considered in this paper the application of 
this approximation is almost trivial. 



2 Application to the Sinai billiard 

We will consider the version of the Sinai billiard jll] consisting of a unit square with a scattering 
disk, having radius R; < 2R < 1, centered on its midpoint. This is a system with infinite 
horizon, for any finite radius R there is a finite number of corridors, though which a particle may 
go without ever bounce onto a disk (these never-touching trajectories constitute a set of measure 
zero in phase space). 

We let the point particle have unit velocity. The trajectory of a particle in the Sinai billiard 
consists of laminar intervals, (bouncing between the straight sections) interrupted by scatterings 
off the central disk. The variable A introduced above is simply the length of the trajectory 
segment between two disk bounces. 



2.1 The distribution of reccurence times 

We introduce the disk as surface of section with phase space coordinates x s . We define A s (x s ) 
as the traveling distance to the next disk bounce. The distribution p(A) is the distribution of 
recurrence times to the surface of section 

p(A) = y J dx s S(A-A s (x s )) , (1) 

where V s = J dx s = AttR using Birkhoff coordinates. 

The relevant time scale is set by the average < A >= J Ap(A)dA which can be computed 
exactly as follows. When we integrate over phase space (or rather the energy surface) it is 
convenient to split up the phase space element according to dx = dx s dr where r is is the coordinate 
along the the trajectory. The total phase space volume may now be expressed as 

r r r A >( x >) f 

V = I dx = I dx s I dr = A s (x s )dx s = V s < A >= AirR < A > . (2) 



On the other hand we have 



V = J dxdy dv x dv v 5{\ - ^vl+v%) = (1 - ttR 2 )2tt , (3) 



leading to the following expression for the expectation value of A 



1 kR 

<A> =2i?~-- (4) 



In ref. [12] we derived the following expression for p(A) in the the small R limit 



12 R 



£<1 



p fl ^o(A)-<j ^ (2e + e( 4_3 01og(e)+4(e _ 1) 2 log(c _ 1) _ (2 _ c) 2 log | 2 _ e | ) e>1 , 

(5) 
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where £ = A/ < A >. Some smearing in £ is required to make the limit well defined fl2]| . 

For finite R there is of course no such simple formula, but the tail (A 3>< A >) exhibit the 
asymptotic law 

KA) ~ ^ , 

where o~(R) is the sum 

a(R) = y £(2qR+^ M -2) . (7) 

where 

S = {q = {n x ,n y )\n y > 0; n x > 0; {n x ,n y ) = 1; q = ^n 2 x + n 2 y < 1/2R} (8) 
is the set of corridors. 

The leading asymptotic behaviour of the function <j{R) is in r cf. found to be a(R) ~ 
l/(47ri? 2 ) and in this limit we get 

2.2 Correlation functions in the BER approixation 

We will restrict ourselves to observables A that are changed only by bounces on the disk. The 
autocorrelation function is e.g. obtained as the time average 

C AA (t) = (A(t + t)A(t )) ta -(A) 2 . (10) 

The BER approximation assumes that there is no correlations if t and t + t belongs to differ- 
ent intervals, that is (A(to + t)A(t )) — (A 2 ) if there is no bounce between t and t + to and 
(A(t + t)A(t )) = (A) 2 otherwise. 

Next let Po(t) denote the probability that the trajectory has not hit the disk between t and 
t + to- We can now write the correlation function in terms of conditional probabilities 

C AA {t) = Po(t)(A 2 ) + (1 - Po(t))(A) 2 - (A) 2 = P (t)V(A) , (11) 

where V(A) is the variance of A. 

The function Po(t) may be expressed in terms of p(A) 



-| />OG />OG 



(12) 



The 1/A 3 decay of p(A) thus implies a 1/t decay of the correlation function. From eq (^) we 
compute the tail of Po (t) 

m ~ _^LI . (13) 

< A > 7T t 

In the R -> limit we get from eqs. @, (§) and (@) 

p o(t)~4p7 ■ ( 14 ) 

TT K t 

For sufficiently small R one can use eq. (0) to obtain an expression for the entire -Po(i) (the 
resulting expression is not particularly nice so we do not display it). For other radii R one can 
use a numerically obtained p(A) to compute Po(t). 

First we will study the observable A = \v x \. Since (|ur| 2 ) = 1/2 and (|« x |) = 2/-7T we get 

C AA {t)=P {t){\-\) , (15) 
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2.3 Correlations in the Lorentz gas 

A closely related problem is the decay of correlations in the associated Lorentz gas obtained by 
unfolding the bounded billiard into an infinite lattice of disks, cf. fig 1. If we kept considering 
the observable A = \v x \ we would just recover the previous result. Let us instead consider the 
observable B — v x . The difference is that the bounces on disks might now change the sign of B. 
The expression ( pi] ) may again be used, with the difference that (B) = so we get 

C BB (t) = Po(t) ■ ({B 2 ) - (B) 2 ) = P (t)/2 . (16) 

It is well-known that this correlation function is related to the diffusion coefficient 

D = lim 1 < (x(t) - x(0)) 2 > , (17) 

i^oo ft 

via the Einstein-Green-Kubo formula |l^] 

D = lim \ [ < v(0) ■ v(s) > ds = lim 2 / < v x (0) ■ vjs) > ds (18) 
t - >0 ° / Jo t ^°° Jo 

where we have made use of the symmetry of our system < v x (0) • v x (s) >=< v y (0) ■ v y (s) >. 
f is the number of degrees of freedom. 

Inserting eq. ( |l6| ) into this formula we obtain the diverging diffusion constant 



U(^-^Llogt, (19) 

< A > 7T 



in agreement with refs. |T^, 



2.4 Comparison with numerical data 

The numerical experiments are carried out through Monte Carlo phase integration (based on a 
subtractive random generator), using 10 7 initial conditions. Data become noisy when correlations 
decrease to the same order of magnitude as the integration error (which is inversely proportional 
to the square root of the number of initial conditions) . It has been pointed out (see for example 
Jl4| and Q) that in principle one might worry about errors due to exponential propagation of 
errors (due to positivity of the Lyapunov exponent and finite precision arithmetics): though 
our time sequences extend beyond the limit suggested by former arguments, the onset of power- 
law (asymptotic) decay is well within the allowed time scale (where the data are insensitive to 
statistical errors, which in turn dominate over the afore mentioned errors). We also remark that 
for hyperbolic billiards Lyapunov-induced errors did not seem to lead to detectable effects, see 

Hi- 

In fig. 2 and 3 we show typical results of our numerical experiments: the numerical correlation 
function indeed exhibits a lit decay for long times. In the same figure we also plot theoretical 
curves obtained from eq. (|15|). The good agreement, even along the preasymptotic region, is 
suggestive of the applicability of BER approximation for such a system. From fig. 2 and 3 one 
can observe that BER expressions, like ( |l5| ) are not able to reproduce preasymptotic oscillations 
in the correlation functions, while being quite close to their envelope, even for moderate time. 
Such oscillations are smaller and smaller in the R — > limit. We comment on the deviation for 
large t in fig 3 in section 3. 

In the small radius limit, eq. (||) suggests a scaling form of the correlation function CAA.nit) = 
g(Rt) where g is independent of the value of R: in fig. 4 we show how the scaling is approximately 
correct for small times, and indeed is well confirmed in the asymptotic region. 

Finally a comparison between correlation functions for A and B observables is presented in 
fig. 5. 
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3 Concluding remarks 

There are two obvious ways of refining the present computations without abandoning the basic 
idea of the BER approximation. 

There is a tacit assumption about isotropy behind the factorization (1 — Po(t))(A) 2 in eq. 
([Ti"l). It is true that the average of \v x \ is exactly = 2/ir but only provided the average is taken 
over the entire phase space. But suppose we restrict the phase space integral to points lying on 
trajectories penetrating deep into the corridors then it is no longer exactly true that {\v x \) — 2/n 
(except in the small R limit). This will lead to a small correction of the prefactor c(R) of the 
asymptotic limit C(t) ~ c(R)/t for observable A but not for observable B. Preliminary results 
indicate that this effect account for the deviation in the large t limit in fig. 3 and fig. 5. 

Secondly, if we have a very long segment into one corridor, than there is an enhanced proba- 
bility that the next segment goes into the same corridor || . Consecutives segments in the same 
corridor should be considered as belonging to the same interval as they are associated with nearly 
the same value of the observable. This effect will provide a small correction to p(A). We will 
return to these problems in future work. 

PD was supported by the Swedish Natural Science Research Council (NFR) under contract 
no. F-FU 06420-303. RA was partially supported by INFN, sezione di Milano and INFM, Unit a' 
di Milano. 
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Figure captions 

Figure 1: The Lorentz gas, obtained by unfolding the Sinai billiard over the plane. 

Figure 2: Experimental correlation function (full line), for R = 0.318. The dashed dotted line 
represents the BER approximation (p"5|), using a numerical p( A). 

Figure 3: The full line represents the numerical correlation function for R = 0.106, the dotted 
line is obtained from (|l5|), using the theoretical small R limit of p(A) (||), while the dashed line 
has been obtained by using the asymptotic expression (p"3|). 

Figure 4: Correlation functions for different R values, plotted as functions of Rt. The dashed line 
corresponds to R = 0.053, the dashed dotted line to R = 0.106 and the full line to R — 0.212. 

Figure 5: Numerical correlation function for observable B (upper full line) and A (lower full line) 
for the Lorentz gas fR=0.159). The dashed lines refer to the asymptotic expression ( |l3| ) fed into 
equations (|l5|) and (^6|). 
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